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Abstract. The two-body Dirac equation with general local potential is reduced to the pair 
of ordinary second-order differential equations for radial components of a wave function. 
The class of linear + Coulomb potentials with complicated spin-angular structure is found, 
for which the equation is exactly solvable. On this ground a relativistic potential model 
of light mesons is constructed and the mass spectrum is calculated. It is compared with 
experimental data. 
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1 Introduction 

It is well known that light meson spectra are structured into Regge trajectories which are 
approximately linear and degenerated due to a weak dependence of masses of resonances on 
their spin state. These features are reproduced well within the exactly solvable simple relativistic 
oscillator model (SROM) [U [21 E] describing two Klein-Gordon particles harmonically bound. 

Actually, mesons consist of quarks, i.e., spin-1/2 constituents. Thus potential models based on 
the Dirac equation are expected to be more appropriate for mesons spectroscopy. In recent years 
the two body Dirac equations (2BDE) in different formulations^ with various confining potentials 
are used as a relativistic constituent quark models PSHEHIBIII21I2QIEIIIZ21I231I231I251I25]- Some 
models are universal \27\ I28j. i.e., describe several states of heavy as well as light mesons. The 
solution of 2BDE is usually obtained by means of perturbative or numerical methods, even in the 
case of simple potentials. As far as exactly solvable models are concerned, only few examples are 
known in the literature which represent versions of two-body Dirac oscillator |29[ 1301 I3T1 132|. [33] . 
Similarly to SROM, they all exhibit Regge trajectories which are linear exactly or asymptotically 
(for large values of angular momentum). But energy levels are spin-dependent, and a type of 
degeneracy is different from that of SROM and experimental data. 

In the present paper we attempt to construct, on the base of 2BDE, the consistent solvable 
potential model which describes as well as possible actual spectra of light mesons. For this 
purpose we start from 2BDE with general two-fermion local potential proposed by Nikitin and 
Fushchich [34]. Following the scheme [351 [28], we perform the radial reduction of general 
2BDE and get a set of 8 first-order ODEs, then 4 of first-order ODEs and finally a pair of 
second-order ODEs. This chain of transformations permit us to write down the set of equations 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
http://www.emis.de/journals/SIGMA/symmetry2007. htmll 

i Here we mention three formulations of 2BDE approach. One of them [U [5] |SJ [7] is built as a generalization 
of the Breit equation [8] [9], two other, [T^[TT][T2] and [13] [14] [15] , originate from Dirac constraint theory. 
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in a compact form and then to impose constraints for general potential which considerably 
simplify the equations. 

In such a way we find a family of exactly solvable models with linear potentials, which 
includes two known examples of the Dirac oscillators [291 [311 [321 [33] and new ones. Then we 
construct integrable extensions of these models describing linear plus Coulomb-like interaction 
of constituents and giving asymptotically linear Regge trajectories. Free parameters can be used 
to provide a spin-independent degeneracy for asymptotics of trajectories. It is remarkable that 
one of the models reproduces exactly the spectrum of SROM. It describes two equal-mass Dirac 
particles bound by rather non-trivial potential. 

Calculated spectra are compared to experimental data for the family of (7r-p)-mesons. 

The Nikitin-Fushchich anzatz as well as our models are rotary invariant but not Poincare- 
invariant, as expected of consistent description of relativistic systems. Thus we propose in 
Appendix the manifestly covariant equations reduction of which into the center-of-mass reference 
frame restores the former general equation with the residual rotary symmetry. In this form our 
models are convenient to compare (and some of them are found equivalent) to other covariant 
two-body Dirac models known in a literature. 



2 General structure of a potential 

We start with the standard three-dimensional formulation of the two-body Dirac equation. It 
has the form: 

{h 1 (p) + h 2 (-p) + U{r)-E}$(r) = 0, (1) 

where #(r) is 16-component wave function (Dirac 4x4-bispinor) of relative position vector r, 

h a (p) = at a -p + mafla = -ia a ■ V + m a [3 a , a = 1,2, (2) 

are Dirac Hamiltonians of free fermions of mass m a , p = — iV, and U(r) is an interaction 
potential. If $(r) is presented in 4 x 4-matrix representation, the operators a a and (3 a act as 
follows: cx.\ <P = a<I>, ct2 $ = <£o! T etc, where a and (3 are Dirac matrices. 

The potential U(r) is a multiplication operator in the position representation, which is in- 
variant under spational rotations and inversion. The most general form of such a potential is 
found in [34] 

48 

U(r) = ^ U A (r)T A . (3) 

A=l 

It is parameterized by 48 arbitrary (complex, in general) functions U A (r) of r = \r\ (we will 
refer to them as partial potentials), and matrices T A are built in terms of Dirac matrices and 
unit vector n = r/r. 

We require of the potential U to be Hermitian with respect to the inner product 

(!?|$) = f d 3 rTr($\r)$(r)). (4) 

Consequently, the equation ([1]) becomes Hamiltonian. This requirement reduces an arbitrariness 
in the general potential ([3]) to 48 real partial potentials. 

Let us construct a Hermitian basis for matrices V involved in the potential ([3]). For this 
purpose we start from the relation for Dirac matrices a = 7 5 <r, where components <Tj (i = 1, 2, 3) 
of the vector a are understood as either Pauli 2x2-matrices or block-diagonal 4x4-matrices 
diag(<7j, Oi). By this a is split into two factors, one of which acts on "particle-antiparticle" 
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degrees of freedom while another one - on spin degrees of freedom. We will refer to them as 
Dirac and spin factors respectively. Similar splitting can be done for arbitrary Dirac operator. 
Now, taking into account that /? is scalar, 7 s is pseudo-scalar, n is vector and a is pseudo- vector 
with respect to 0(3) group, we write down the basis for T-matrices in the product form 



r e 

To 



{I, 01, X S, 

{/,/3 1 ,i/3 2 ,i/3 1 /3 2 }x 7 |xT, 



7i 5 7l x r c 
r 1 

J- eo| i<_>2 ' 



(5) 
(6) 



where I is the unit operator, S and T are sets of scalar and pseudo-scalar spin factors 

S = {S {1) , a = l,2,3} = {I,<7i-<7 2 ,(<7i-n)(o- 2 -n)} (7) 



T = {T (i 



1, 2, 3} = {a 1 ■ n. rr 2 ■ n, (n, tr 1 , cr 2 )}. 



(8) 



The matrices T in ([5]), © are grouped as the even-even operators r ee , odd-odd r oo , even- 
odd r co and odd-even operators r oc , following the Chraplyvy classification [36]. All other 
0(3)-invariant Hermitian matrices built with Dirac matrices and n can be expressed in terms 
of ©, ©• Correspondingly, the potential (j3J) can be split into four terms U ec , U 00 , U co , U oe . 

It is convenient to use block-vector representation for the wave function instead of matrix 
one as follows 



<P(r) 



(r) 
(r) 



(r) 



*++(r)' 

*+-(r) 
<P__(r) 



here #_|_ + (r) is 2x2-spinor matrix which represents large-large component of wave function 

while (r), # |_(r) and # (r) are large-small, small-large and small-small components 

of $(r). In this representation the 48 of partial potentials [/^(r) in © are collected in Dirac 
multiplets with common spin factors as follows 



i=i 



A ee + A 00 , = fl co + ft oc (here we omit a superscript i) and 



u 22 
n 1783 

U c/44 







^23 



1^23 



Ti- 



ll 



wi 2 

Wl2 









w 34 

W34 



ft, 





n Wis 

U Ty 24 

Wis 

W" 2 4 U 



(9) 



(10) 



(11) 



where matrix entries U a p are real while W a p = X a p + iY a p are complex functions of r (the star 
" * " denotes a complex conjugation). 



3 Radial reduction of 2BDE 

Following [351 171 [28] we choose #(r) to be an eigenfunction of the square j 2 and the component js 
of the total angular momentum j = rxp + s = — ir x V + \{tri + <r 2 ) and of the parity P. In 
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a block-matrix representation we have 



#(r) = - 
r 



0\ 



\<t)i{r)<3f A {n) + i(/> 2 (r)^°(n) ' 
3 (r)^-(n) + 4 (r)^+(n) 
5 (r)^-(n) + 6 (r)^+(n) 
i0 7 (r)^ A (n) + i0 8 (r)^°(n) _ 

-+(^~,^+) forP=(-) J ' 



for P= (-) j±1 , 



Here ^ A (n) is an abbreviation of the 2x2 matrix bispinor harmonics ^^{n) (in conventional 
notation |37t [29l I27j ) corresponding to a singlet state with a total spin s = and an orbital 
momentum £ = j; harmonics & (n), ^ _ (n), <3f + (n) correspond to triplet with s = 1 and 
I = j,j + l,j — 1. Then for j > the eigenstate problem ([1]) reduces to the set of eight first- 
order differential equations with the functions (f>i(r), . . . , 4>%{r) and the energy E to be found. 
Using properties of bispinor harmonics (their explicit form is not important here; see |35|, [28] 
for one) 



J dnTr(&? &k) = 6 ik> 



j 2 W = j(j + l)W, i = 0,1,..., 
j 3 W = H&, n = -j,...,j, 
£ 2 & = £(£ + l)W, £ = j,j±l, 
s 2 <3f = s(s + l)W, s = 0,l, 

and defining a 8-dimensional vector-function <t>(r) = {<pi(r),. 
the matrix form 



H(j)^.+V(r,E,j) 



Here 8x8 matrices H(j) and 

V(r,E,i) = G(i)/r + m + U(r,i)-^l 



, 08(r)} we present this set in 

(12) 



(13) 



possess the properties: H £ Re, H T = — H; = V; I is 8x8 unity; the diagonal matrix m and 
matrices H(j), G(j) are constant (i.e. independent of r), and U(r, j) represents the potential ([3]). 
The operator in the l.h.s. of (|12|) is Hermitian with respect to the inner product: 



(v|/|d>) 



dr (v|/t( r )d)(r)) 



induced by ([!]). In the case j = components 4>2 = </>4 = </>6 = <5 
the problem ()12|) reduces from 8 to 4. 

Uniting bispinor harmonics in two doublets of opposite parity 



'8 = so that a dimension of 









^0 ^ 


e = 
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with the properties 

CT | • (T2 6 = T 6, 

<t\ ■ no = R e, 
a 2 no = — o"3 R T e, 
(<7i • n)(er 2 • n)6 = -cr 3 e, 
(n, (Ti, <r 2 ) 6 = — 2icTtR T e, 



<t i • cr> e = e, 

(Ti ■ ne = R 6, 

<T2 • n e = — Ro"3 6, 

(<ri • n)(er 2 • n) e = -Rf7 3 R T 6, 

(n, <ti, <x 2 ) e = 2i J Rcr-f e, 



(14) 



where 2x2 matrices o>, en, r and R are defined as follows 

ct t = i(l + cr 3 ) 
_3 n 

R 



1 





1 



-3 
1 

JTT 

2j + r 



A B 
-B A 



B 



j 



a i = 10 - °"3) = 

(R G 0(2)), 

c = vTcJ+i), 



2j + r 

we present all the matrices involved in the equation (112f) in a block-matrix form 

m = diag(m + l, m_l, — m_l, — m+l) 
(here m-t = mi ± m 2 , and I is a 2 x 2 unity) , 





(i+a T )R T 

C?0"3+cr T )R 





R T (j£7 3 +fT T ) 







R(j+a T ) 



H 





a 3 R T 
R 




-Ra 3 



-R T 







H 



R(j-h7 T )cr 3 

R T 0'o-3+o--r) 
o 

o -0-3 R T 
Ra 3 
R T 

-R 



R(j+<7 T ) 



R t (j^ 3 +^t) 

— R T 

R 

0-3 R T 
-Rcr 3 

(i+(T T )R T (j(J 3 +CT T )R 





(ja 3 +tT T )R (j+<r T )R T 

-R 
R T 

Rcr 3 





(io" 3 +cr T )R 

(i+a T )R T 




for P = (-) 







7±1 







R T (jo-3+o-t) 
R(j+a T ) 





for P = (-) J 



In order to present the interaction potential (|9j)- (llip in the block-matrix form we have to 
calculate in this representation the form of spin operators ([7]), (jHJ) using the relations (I14p . One 
obtains 



S = giag(S, E, E, S), T = giag(T, -T+, -T+, T) 
S = giag(E, S, S, E), T = giag(T+, -T, -T, Tt) 
with the following 2x2 matrix blocks S, E, T 



for P 
for P 



i 


1 


2 


3 




J 


<T\ ■ (72 


(tri • n)(<r 2 • n) 


S (i ) 


1 


T 


-ct 3 




1 


1 


— R<73 R T 



i 


1 


2 


3 




rr | • n 


cr i ■ n 


(n,a- 1 ,a 2 ) 




R 


-Ra 3 


-2iRo- T 
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In the case j = all 2x2 blocks in matrices of the equation (|12p must be replaced by their 
left upper entries calculated at j = 0. 



4 Reduction of 2BDE to a set of second-order ODEs 



It turns out that rankH = 4 (2 for j = 0). In other words, only four equations of the set (|12p 
are differential while remaining ones are algebraic. They can be split by means of an orthogonal 
(i.e., from the group 0(8) ) transformation 

OHO T etc, 



cD = 0<D, 



Oi 



Oi 



I 



H : 

<r 3 R T 









a 3 R T 





-Ft 



-cr 3 



where O 



2 = ^ 







V2 



1 

71 



2 Oi, 






1 





-1 " 


-1 


1 





1 


1 





1 





1 


1 


-1 





1 





-1 


-1 





-1 


1 


1 






7±1 



for P = (-) 



for P = 



This transformation reduces the matrix H to the canonical form 
H = 2 



JO" 




o r 





J = 


-10 



(15) 



where J is a symplectic (nondegenerate) 4x4 matrix. Other items of (I12p - (|13p take the form 

m± 



U 






(7 | 




" 

















Coi 



























Un 


W12 


Wig 


W14 


Wl2 


u 22 


w 23 


w 24 


Wl3 


W 23 


U33 


W34 












Wl4 


W24 


W34 


U44 



m 







m± 



-m- 







forP = T(-) j , (16) 



(17) 



here U Q/ g are real and \N a p are complex diagonal 2x2 blocks (related linearly to the former partial 
potentials U a p, W a p in (fT0]) - (fTT]) ) which are different, in general, for the same j but opposite P. 
Equivalent amount of real partial potentials present in (|17p is equal to 32 for P = (— ) J±1 and 32 
for P = (— y , but the only 48 among them are independent. 

Let us present 8-dimensional vectors and matrices involved in (|12p via 4-dimensional blocks 



0! 

d> 2 



V 



Vn v 12 
V 2 i v 22 



and (I15p for H. In these terms the set (|12p splits into differential and algebraic subsets 

2J*' 1 + V 11 * 1 + \/ 12 * 2 = 0, (18) 

V 21 0! + V 22 2 = 0. (19) 

Eliminating <t> 2 from (I18p by means of ()19p yields the purely differential set for 4-vector Q>i 



0. 



(20) 
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The vector <t>2 is then determined by the algebraic relation <t>2 = — AV2i4>i. Here 

A = [V 22 ] ~ l , V x = (V n - V 12 AV 21 )/2. (21) 
Next, we present the 4-vector 0i in a 2+2 block form 



*i(r) 



0! 
02 



Vn v 12 
V 2 i v 22 



An A12 
A21 A22 



(22) 



then eliminate 02 and arrive at the set of second-order differential equations for 2-vector 0] 



d 
dr 



V21 + Vn \ 0i = 0. 



(23) 



It follows from (Q2J), CE]), ([T7D, (JH]), ([22]) that the matrix V 22 is diagonal. Thus one can 
perform the transformation 

V = 0i/7V^, L(E) = y^l(E)^ 2 , 
reducing the equation (|23p to a normal form 



where 



l(E)MJ 

3 

F : 

Z 



d 
dr 



+ F 



d 
dr 



FT + Z \ V = 



(24) 



'V 22 V 



22 V12 



V 



22 



and 



/T7- 1 ^ V 22 V 22,.A 1V 22 3 /V^ 2 

V 22 Vii VV 22 + - I + r#F T ) + 



V22 V22 



2 V 



22 



4 V V 



'22 



(25) 



If F is Hermitian then the first-order derivative term is absent in the operator \-(E) and the 
equation (|24|) takes the matrix two-term form 



A2 

_ + Qj>vl/ = with Q = Z-F 2 -F' 
dr z 



which is useful for a search of solvable examples of 2BDE. 



(26) 



5 Solvable oscillator-like models 

A construction of the matrix V 22 d22]), §Z\§ involved in the equations ([20]), ([23]), ([M]) and ([26]) 
is a source of non-physical singularities which very complicate an analysis of these equations. 
Moreover, singularities may make a standard boundary value problem (with constraints at r = 
0, 00) incorrect. We note that in the free-particle case V22 is a constant (free of r). Let us require 
of V22 to be so in the presence of interaction. Sufficient conditions for this are the equalities 



U 



22 



U 



33 



w 



23 



w 



24 



W 



34 







(27) 



for entries of the matrix fjlTf) leading to a rather general potential of the form ([3]) with 18 arbitrary 
real partial potentials. Four of them compose the even-even part U ee of the interaction, four 
other constitute the odd-odd U OCn and the remaining 10 ones form potentials of the U oe and U eo 
type. Due to the constraints (|27|) both the £7 CC and U QO interaction terms contain the factor 
1 — cri • cr 2 . It cancels the interaction on a triplet part of the wave function $ and leads to a 
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dubious bound state spectrum of poor physical meaning. Thus we discard these potentials by 
means of the conditions 



U 



ii 



U 



11 



w 



14 







for every of the P = —(—)•' and P = (— )•? parity cases. 
The following requirement 



ReW 



12 



ReW 



13 







(28) 



causes the equality Im F = for (|25p and thus = F (since F is diagonal, due to the equali- 
ties (|27p ). By this we provide the two-term form (|26p for the wave equation. On the other 
hand ten partial potentials contained in the entries W12 and W13 of the matrix (|17j) reduce, due 
to (|28p . to 6 ones which form interaction terms of U oe and U co type. Among them the term 
Uo(r)(ai — a 2 ) ■ n with an arbitrary function Uo(r) is purely gauge and can be compensated by 
the phase factor exp {—if drUo(r)} of W. The resulting anzatz for a potential is as follows 

U = i{U 1 (r) + C/ 2 (r)7i 5 7 2 5 } PMai-cx 2 ) ■ n + U 3 (r) ( 7l 5 + 7 |) fa x <x 2 ) ■ n 
+ i{U 4 (r)Paf - U 5 (r)f3 2 ^} (°"i ~ °"2 + i"i x a 2 ) ■ n, 

where U\(r), . . . , Us(r) are arbitrary real functions. 

After aforementioned simplifications the matrices Q+ and Q_ involved in the equation (|26p 
and referred to the P = — (— ) 3 and P = (— ) J parity cases respectively have the form 



Q 



1 



± 



E 



rn 



+ 



E 



E 



/f± + 



d 
dr 



m_ 
~E 

/t± 
c 



~E 



2/ T ± + 



dr 



m 



-f(h±±fl)-u± 



u± + u± 
C 2 



a 1 



fl 9") 



where five functions 

h± = T(Ui + U 2 ) - 2U 3 , fi = U 2 - U h 



u± = -2(U A ± U 5 



represent equivalently partial potentials. They are a convenient choice when looking for solvable 
models. 

Models I. The choice u = and f 2 = a 2 r 2 , where a = const, leads to four possibilities 



la) fu 
lb) f U 



-/t- = -h 

-/t- = h = 



—ar, 
-ar, 



Ic) /t+ = /?- 
!d) /t+ = /t- 



which correspond to the following sets of partial potentials (only non-zero ones are shown) 
la) U 2 = ar, Ic) E/i = -U 2 = U*- 1 - 



'2 
-U 2 



lb) U! = ar, Id) Ui 
All cases lead to the oscillator-like equation 



3 - 2 
^3 = 



ar, 



7 ar. 
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where D± is a constant 2x2 matrix. In the case la) this matrix is diagonal, 

la) D±= diag{,5 T± ,o i± }, (30) 
*T± = ±3, *i± = =Fl forP = T (-) i - (31) 
The energy spectrum can be found from the algebraic relation 
1 / Trii_\ ( m 2 



E-^UE- — )=\a\(2n + 3)-a8 l± , n = j + 2n r , (32) 



where n r = 0, 1, . . . is a radial quantum number. Given j and n r we have 4 equations determining 
4 positive values of energy E^± 



a=l 



In case lb) the matrix D± is not diagonal, 



3 -2m T C/E 
-2m T C/E 1 



forP = T{-) j - (34) 



It can be reduced, by means of 0(2)-transform, to the form (|30p with eigenvalues 

lb) o T± = 2 + ^l + (m T C/E) 2 , 5 l± = 2 - y/l + (m T C/E) 2 for P = =f(-) 3 '. (35) 

The set (j29[) splits into two oscillator-like equations. Since S's (|35p are energy-dependent, the 
spectral conditions (|32p appear in this case as irrational equations which can be reduced to 
a fourth-order algebraic equations (third-order if mi = 7712) with respect to E. 

In cases Ic) and Id) the matrix D± is equal to either (f30|) - (f3"Tj) or (j34"|) for P = (— case and 
to another (complementary) one for opposite parity. The corresponding eigenvalues 



Ic) 5 t+ = 2 ± ^1 + (m_C/ J B) 2 , 5j_ = 1 ± 2, 
Id) 5 t+ = 1 ± 2, Jj_ = 2 ± v^l + (m + C/E) 2 

are to be used in the spectral formulae (|32|) or (f33l) . 

If mi = 7?i2 = 0, the constants 5^± do not depend on E and j and the equation (|32p simplifies 
to the explicit formula for the energy spectrum 

J E I 2 ± = 4[|a|(2n + 3)-a5 I± ]. 

In this case energy levels E^± build up in the (E 2 , j)-plane into parallel straight lines (so called 
Regge trajectories) with the slope rate 8|a|. All the trajectory can be labeled unambiguously 
by the triplet of quantum numbers n r , |, ± but some trajectories corresponding to different 
numbers may coincide (i.e., degenerate). Examples of Regge trajectories for massless models 
Ia)-Id) are shown in the Fig. [TJ 

Models II. The choice / T = and ff = u 2 = a 2 r 2 leads to four possibilities 

Ha) fi = ±u± = —ar, lie) fi = u± = —ar, 
Hb) fi = =fu± = -or, lid) fi = -u± = -ar, 

which correspond to the following sets of partial potentials (only non-zero ones are shown) 

Ha) Ux = U 4 = \ar, He) U t = U 5 = \ar, 

lib) U x = -U 4 = -\ar, lid) Ui = -U 5 = -\ar. 
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la). 





Figure 1. Rcgge trajectories from massless 
models I and II. 




o - 



II). 




-ot+ 

■ 01+ 



-H,t+ 
-n |+ 



• Of- etc., where n =1,2, ... 

■01- 



This choice leads to four solvable models described by the oscillator-like equation (|29p with 
some matrices D. Here we consider the massless case only. All models IIa)-IId) have the same 
spectrum 



a > 



E 



1 + 



: a < : 

4a(3j + 4n r + 3), E* + = 4|a|(j + 4n r + 3), 

£ 2 + = 4a(j + 4n r + 2) , £ 2 + = 4|a| (3j + An r + 4) , 

= 4a(i + 4n r + 3) , £f_ = 4|a| (3j + 4n r + 3) , 

E\_ = 4a(3j + 4n r + 4), E\_ = 4\a\{j + 4n r + 2). 



In contrast to models I, here the Regge trajectories are not parallel: two families are sloped 
by 4|a|, two other - by 12|a| (see Fig. [T|). 

Models III. The conditions / T - = u+ = 0, / T 2 + = ff = u 2 _ 
one example of which corresponds to the potentials 



a 2 r 2 yield four solvable models 



far, 



Uo = Uz 



U, = U 4 



r ar. 



These models (in the massless case) lead to sequences of Regge trajectories with three different 
slopes: 4|a|, 8|a| and 12|a|. 

Two of the above models are known in the literature. Model la) is equivalent to the Sazd- 
jian pseudo-scalar confinement model [29]; it coincides with one of versions of two-body Dirac 
oscillator proposed in [32]. Supersymmetric aspect of these models has been studied in [30| [33] . 
Model lb) is a generalization of another version of Dirac oscillator [31] to the case of different 
rest masses. 

While having no direct physical meaning, the energy eigenvalues of the models possess some 
important features of actual meson spectra which we discuss in the next section. 
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6 Light meson spectra 

Characteristic features of mass spectra of light mesons (consisting of u, d and s quarks) can be 
summarized roughly in the following idealized picture [381 139"! [71 [28] : 

1. Meson states are clustered in the family of straight lines in the (E 2 , j)-plane known as 
Regge trajectories. 

2. Regge trajectories are parallel; slope parameter u> is an universal quantity, u = 1.15 GeV 2 . 

3. Mesons can be classified non-relativistically, as (n 2s+l £j\ -states of quark-antiquark system 
(n = n r + £ + 1 where £ and n r are the orbital and radial quantum number) . 

4. Spectrum is .^-degenerated, i.e., masses are distinguished by £ (not by j or s) and n r . 

5. States of different £ and n r possess an accidental degeneracy which causes a tower structure 
of the spectrum. 

It follows from items 1-3 that there exist 4 principal (n r = 0) trajectories, one of which 
includes singlet states (s = 0, j = £), and three others collect different triplet states (s = 1, j = 
£,£± 1). Each of principal trajectories heads a sequence of daughter trajectories (n r = 1, 2, . . . ). 
Item 4 means that in the (E 2 ,£) -plane the four parents degenerate completely (this is concerned 
also with daughters of the same n r ). Then energy levels of q-q system can be described by 
a formula 

E 2 ^Lj(£ + Kn r + (), (36) 

where a constant £ depends on a flavor content of mesons (C ~ 1/2 for (7r-p)-family of mesons; 
it grows together with quark masses). Finally, item 5 restricts the constant k to an integer 
or rational number [39]. One frequently puts k = 2 (oscillator-like degeneracy) [TJ [2j [31 HQ], 
rarely k = 1 (Coulomb-like one) [4RI42j. Fig. [2] illustrates the case of (£+2n r )-degeneracy where 
corresponding Regge trajectories are plotted in planes (E 2 ,j) and (E 2 ,£). 

Actual spectra of light mesons differ from the idealized spectrum (I36h . A finite number of 
mesons are known, some meson states are ambiguously identified by quantum numbers, Regge 
trajectories are not quite straight (especially, in their bottom) [43] degeneracy is approximated 
(~ 5 -r 6% of oS) etc. Thus the equation (|36p represents light meson spectra as a rather rough 
approximation, and its subsequent correction is implied. 

Most of features described above are characteristic of spectra of the models Ia)-Id) (see 
Section 5). Indeed, in the massless case Regge trajectories are straight and parallel. The 
trajectories |+ consisting of energy levels E^ + of parity P = — (— ) 3 can be treated as a singlet 
trajectories corresponding to s = 0, i = j (and labeled by "A"; see Section 3). Similarly, triplet 
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quantum numbers s = 1, I = j,j =F 1 can be prescribed to the trajectories |+, |— (labeled 
by "0", " = F"). In these terms the spectrum of the equal-mass (mi = mi = m) model la) is as 
follows: 

E 2 = 8\a\(£ + 2n r + 3/2) + 8a(2s - 3/2) + Am 2 
which agrees with results presented in [291 [30], [32| [33] . For the model Ic) we have 
E 2 = 8|a| (£ + 2n r + 3/2) + 8a(s - 3/2) + 4m 2 

(the only difference consists in a coefficient at s). 

All models Ia)-Id) reveal degenerated trajectories. But the degeneracy is not a is one. In 
order to provide £s-degeneracy one should have a model possessing additional (3(3)-symmetry 
with corresponding conserved generators [44]. The total spin s = \((T\ + <r 2 ) does not match 
for this role since even the free two-body Dirac Hamiltonian does not commute with s. Below 
we propose generalizations of the I-type models which reveals an approximated £s-degeneracy 
at some values of free parameters. 

Model IV is an integrable extension of the model lb). We put u = and 



/t± = Tar + (x± - l)/r, = -ar + x/r, 



(37) 



where x+, x_ and x are arbitrary constants. So the model includes, together with a slope 
parameter a and rest particle masses mi, 7712, six arbitrary parameters. 
The choice (|37p corresponds to the following partial potentials 



x 4 



ar 



x_ + 2 X 



U 2 



2x 



Ar Ar 
The matrix Q± in the wave equation (|26p has the form 



4r 



(38) 



Q 



± 



^ 1 



~E 



2„2 



c 2 + x(x±i) 



+ 2 [±a 



2r 2 



-?n T C/£ 



-m T C/E 




(39) 



The matrix in r.h.s. of (J39j) and thus the Q± itself can be diagonalized, similarly to case lb). 
Then the equation f|26|) splits into two oscillator-like equations with j- and i^-dependent free term 
and coefficient at 1/r 2 term. Spectral conditions which in general are cumbersome irrational 
equations cannot be solved explicitly. In a massless case we get exact expressions for Regge 
trajectories 



E^± = 8\a\(k u + 2n r + |) T 8a(x ± + ±), 



k 



Ef± = 8\a\(k l± + 2n r + |) - 8a( X ± +\ 



2)i 



k 



(j + i) 2 + x±(x± 
0' + ^) 2 + x(x ± 1) 



1) 



1 

2' 



1 

2 1 



The Coulomb-like terms in the potentials (|38p bend Regge trajectories (especially in the 

j — >oo 

bottom) leaving their asymptotics (at j — ► 00) rectilinear (since ~ j + 0(1/ j)). When 

choosing x + = x, x_ = — x — 2 (where x is still arbitrary), trajectories become degenerated 
asymptotically in the plane (E 2 ,£). This model satisfies approximately item 4 of the properties 
of meson spectra but leads to top-heavy meson masses. 

Model V is another extension of lb) with the same amount of free parameters as in model IV. 
We put 



u± = 2x±/r, /|± = T{ar - (x T l)A\h fi = ~ar + x/r, 
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which correspond to the potentials 



X 

ar , 

r 



U 2 = 0, 



lh 



l/2_ 

r 



lh 



2r 



lh 



^4 



2r 



The matrix Q± is as follows 



Q 



l 



± 



E 



m j 



E 



E 



m_ 
~E 



a 2 r 2 + a(2x±l) 



c 2 + x (x±i) 



2 ±- 



m- 



E 



XTl/2\ * 



2x± C 
C 



Again, an evident diagonalization of Q± leads to a split pair of oscillator-like equations and 
spectral conditions similar to case IV. In a massless case we have the spectrum 



Ef = 8\a\(k t + 2n r + §) - 8a( X ± \), 



where 



k 



1 



C 2 + (XT W + 2|x|(|x| ± v^ + x 2 ) - 1/2 « j ± H + 0(l/j) 



and parity indexes "±" are implied. Regge trajectories are asymptotically linear but there is no 
a choice of parameters providing the £s-degenegacy. 

Models IV and V represent maximal solvable extensions of model lb) found in this work. 
There exist extensions of other model of family I which reflect several features of light meson 
spectra but are exactly solvable with less number of free parameters than models IV and V. 
Here we consider one example which has close relevance to a meson spectroscopy. 

Model VI is the integrable extension of model Ic). We put u = and 



ft+ = -ar + (x - l)/r, ft- = ft 
which corresponds to the potentials 
x + x ~ I s 



-ar + x/r, 



u 3 



ar 



ar + 



2r 

x + x- 1 

2r 



lh 



1 / , *-3 X -l 
■ — ar H 

2 V 2r 



The matrix Q + is identical to that (|39p of model IV and a treatment is the same. The matrix Q_ 
in the present case is diagonal 



Q = diag{Qf_,Q;_}, 



Qi- 



E 



mfj. 
~E 



E 



m_ 
~E 



oV + o{2(x±l) + l} 



c 2 + x (x + i) 



In the case of equal particle masses, mi = ni2 = m, we get the energy spectrum explicitly 

E 2 + = 8\a\(k n + 2n r + §) -8a(x + §) +4m 2 , fc T+ = ^/ ( j + ±) 2 + x(x- 1) - ±, 
Ef + = 8\a\(k l+ + 2n r + §) - 8a( X + ±) + 4m 2 , = sj(j + ±) 2 + x ( x - 1) - (40) 

£f_ = 8|a|(fcj_ + 2n r + |) - 8a( X + ± ± 1) + 4m 2 , fcj_ = ^(j + ±) 2 + x ( x + 1) - ±. 



Regge trajectories are asymptotically linear and, if x = X , ^s-degenerated (see Fig. [3|). 
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Figure 3. Regge trajectories from model VI; mi = = 0, x = x = 1/2. 

Let us try to describe the spectrum of light mesons of the (7r-p)-family by means of equa- 
tions (|40p using four arbitrary parameters a, x, x and m as adjustable parameters. We note 
that the intersection of the p(770)-trajectory (which is Oj— in our terms) with the £? 2 -axis in the 
(E 2 , j)-plane is negative while it follows from (j40]): E? (j=0, n r =0) > 4m 2 . On the whole one 
can achieve a qualitative and partially numerical agreement of the model with experimental data 
when supposing that the meson masses squared M 2 are related to E 2 as follows: M 2 = E 2 — c 2 
where the parameter c 2 > is common for all stated. Similar picture arises within other 
models, in particular, in [24] where the meaning of the constant c is discussed. Alternatively, 
one supposes m 2 < (instead of the use of c) in some SROM potential models [21 05]. This 
is not acceptable here since imaginary rest masses break the Hermicity of the two-body Dirac 
Hamiltonian ([I]). 

It is shown in Fig. [4] the spectrum of the (vr-p)-family plotted with the data [46] and corre- 
sponding Regge trajectories built with optimal values of parameters: a = 0.145 GeV 2 , x = 0.227, 
k = 0.342, c 2 = 0.346 GeV 2 . Principal trajectories fit majority of meson states. Some radially 
excited levels, however, are situated far from daughter trajectories. 

In the case x = x = trajectories become linear and degenerated exactly, and formulae (|40h 
drastically simplify 

E 2 = 8|a| (£ + 2n r + f ) + 4(m 2 - a). 

Except for a definition of an additive constant this formula describes exactly the spectrum of 
the relativistic system of two spinless particles harmonically bound (SROM) [H [H [3]. In the 





2 The parameter c 2 absorbs in (|40|) the term 4m 2 which then is set to zero. 
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Figure 4. Spectrum of it- and p-mesons and optimal Regge trajectories from model VI. Reliable data 
are inscribed in bold, doubtful - in italic. Thin horizontal lines denote measurement errors. 

present case this simplicity is achieved owing to a rather nontrivial choice of the total potential 

U = l-ar /3i&(7i " tIK^i + <x 2 ) • n + ^(3Ml\ + tD^i - <r 2 ) ■ n 
2 4r 

+ \ ( ar + h) ^ + 1 ^ (Tl x a2) ' n ' 

An exact ^s-degeneracy indicates the existence of additional 0(3)-symmetry which cannot 
be generated by components of the total spin s = \{(T\ + <T2). A search of relevant conserved 
quantities is beyond the scope of this work. 



7 Summary 

In this work we consider a two-body Dirac equation with general local (in the position represen- 
tation) potential found by Nikitin and Fushchich in [34] . It is parameterized by 48 real functions 
of r and is presented here in a matrix form. 

Owing to 0(3)-invariance of 2BDE it is reduced (via a separation of variables) to a set 
of eight first-order ODEs for radial components of the wave function. Then using a chain 
of transformations we eliminate some components of the wave function in favor of other ones 
and arrive at the related pair of second-order ODEs. Coefficients of these equations have, in 
general, poles at some energy-dependent points te which are absent in an original potential. 
These singularities may complicate an analysis and a solution of the problem. 

The structure of second-order reduction of general 2BDE suggests a wide class of potentials 
(parameterized by 14 arbitrary functions) for which the problem is free of non-physical singulari- 
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ties. Within this class a family of exactly solvable models is found which generalize known 
two-particle Dirac oscillators EH EH E21 E3] ■ In particular, two of these (models VI and V) 
are 6-parametric integrable extensions of the oscillator model [31 j . Regge trajectories following 
from these models have parallel rectilinear asymptotes but are curved in their lover segments. 

The 5-parametric model VI is used as a solvable potential model of light mesons. Special 
choice of parameters leads to linear Regge trajectories which possess an exact £s-degeneracy 
and the accidental (j+2ra r )-degeneracy. In other words, the model restores the idealized meson 
spectra generating by SROM [H El E]. The corresponding two-fermion interaction potential 
turns out to be surprisingly intricate. 

Unfortunately, the model fails to describe properly the spectrum of lightest mesons as it 
overestimates meson masses squared M 2 by certain (common for all states) constant c 2 . Up to 
this discrepancy the description of (7r-p)-family is adequate and has been realized. The choice 
of parameters in this case differs slightly from that of the completely degenerated model so that 
the Regge trajectories are somewhat curved in the bottom. The fit is good for most of orbital 
meson excitations and worse - for radial ones due to shortage of daughters following from the 
model. It follows from this fact that the accidental (j+2n r )-degeneracy inherent approximately 
or exactly in this and many other models [H EJ El [40] is less adequate to actual meson data than 
the degeneracy of (j+re r )-type. 

Despite the 2BDE used in this work being 0(3)-invariant and not truly relativistic equation, 
it can be considered as some covariant equation reduced in the center-of-mass frame of reference. 
In Appendix we constructed explicitly the Poincare-invariant equation for 2BDE with potential 
of general form and proved its unambiguity. 



Appendix. Covariant form of 2BDE 

The two body Dirac equation considered in this work is 0(3)-invariant but not Poincare- 
invariant. Below we construct the manifestly covariant equations reduction of which in the 
center-of-mass reference frame restores the 2BDE ([1]) with arbitrary potential of the form ([3|). 

We start with a free-particle system. Following by the constraint formalism |10[ \TT\ fl~2| [T3"l 
[T4l [T5] it is described by the pair of covariant Dirac equations 

(7a -Pa -m a )$ = 0, a = 1,2, (A.l) 

where 7 a • p a = ^/aPafi, the particle 4- momenta p afJ , (a = 1,2; n = 0, . . . , 3) are conjugated to 
particle positions Xa, and the timelike Lorentz metrics \\rju,v\\ = diag(+, — , — , — ) is used. 

Before introducing an interaction it is convenient to represent the equations (jA.ip in a collec- 
tive form. For this purpose we perform the canonical transformation (xa,p afl ) — > (X^, P^^x^jP^) 

x = xi - x 2 , P = P\+P2, M = \fP 2 ~, 

p = pi-e(M)P, X = £x 1 + (1-0x2 + ^(P-x)P, P = P/M, 

where £(M) is an arbitrary function, and introduce the operators 

tt free = 7 Jj(7i • Pi ~ ™i) + 7i'(72 • P2 ~ m 2 ) 

= 7i72 M + 72(7i -p-mi) +7i(-72 -p- m 2 ), 
= \(lx ■ Pi + rai)(7i • Pi ~ m i) ~ 5(72 • P2 + m 2 )(72 • P2 - m 2 ) 
= \{pi-p 2 2-m\+ml)=P-p-iy(M), 

such that pi hcc ,%] = 0; here 7^ = -y a ■ P and v(M) = \ (mf - m| + (1 - £{M))M 2 ). In these 
terms the equations (|A.1|) take the equivalent form "K free <?> = and %<!> = 0. 
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At this point we choose the mixed representation for the wave function: & = $(P, x). Then 
it follows from the constraint K<P = 

$(P,x) = exp U'W-A $(P,x ± ), 

where x" = P ■ x ^ x° and x^ = (6g - P^P v )x u ^ (0,r) thus $(P,x ± ) $(E,r) where 
E = Po|p=o = M\ p=0 . In other words, this constraint suppresses the ^''-dependence of $ and 
eliminates the relative time x° from the center-of-mass description. Equivalently, one can choose 
v(M) = 0, then $ is free of x". We note that the choice of v(M) or £(M) does not affect the 
form of the operator J{ free and so it is gauge-fixing. Thus we do not distinguish the functions $ 
and any more. 

In the center-of-mass (cm.) frame of reference (where P = 0) we have 
!K frcc <2>(P,x ± ) {E - hip) - h 2 (-p)} #(r) = 0, 

where p = p 1 = —p 2 , Pa = la = 7"|p_ , an d operators h a (p) are defined in Up to the 
factor —fiify this equation coincides with 2BDE JT]) in the case of free particles. 

In the general case the interaction term V<P appears in the r.h.s. of the equation J{ free <P = 
(instead of zero) , and we have the following set of equations 

P ■ p V = and :K frcc ^ = V*. 

Since a free-particle operator !K free is Poincare-invariant the potential V must be so. Besides, it 
must obey the equalities 

Vc.m. = V| p=0 = lhfaU(r) (A.2) 

and [V, P • p] =0. These requirements allow us to construct unambiguously the covariant 
operator V for arbitrary potential of the form (jSj). 

It follows from (|A.2|) . ([3]) that the operator V c . m . consists of the sum of partial potentials Ua(t") 
multiplied by matrix coefficients Pip2^A- Both factors in every term of the sum are 0(3)- 
invariant. We construct Lorentz-scalar counterparts for them separately. 

Let p = \J — x\ and thus p c . m . = \r\ = r. Then functions Ua(p) are Lorentz scalars by 
construction and UA(p) cm =Ua{t). 

©-© is a basis for matrices Pifh^A- PiPi^ee = r ee , PiP 2 T OQ = T OQ) PiP 2 T eo = ir co , 
P\P 2 T oc = ir oc (the factor "i" is not important in this regard). First of all it is convenient to 
express er-matrices in the spin factors ([7]), (jHJ) of the basis {[5l) — (EJ) via 7-matrices as follows: 
a a = JaPa^fa- Then every basis element is a product of factors 

I, 7a, Pa, 7l'72i 7 a '«, (»*, 7l,72) 
which are scalar or pseudo-scalar with respect to the group 0(3). We construct the Lorentz 
scalar or pseudo-scalar counterparts of these matrices as follows 

I —* I, la —* 7a, Pa — » 7a, 7l " 7 2 ~ * -ll'l2, 

7 a -n — » -ia-n±, (n, 71,72) — » ~(A «,7i, 72) = -^A^ M «^7^, (A. 3) 

where n = x/p and e^Ax is absolutely antisymmetric pseudo-tensor (£0123 = — !)• 
The resulting covariant potential V has the form 

V(x,P) = Y?* =1 Ua(p)Za, 

where matrices Ea are built with P\P 2 ^ A by the replacement (IA.3j) (with an unchanged order 
of co- factors). 

The unambiguity of this construction is obvious: if there exist two potentials which coincide 
in the cm. reference frame then their difference is zero in any reference frame, due to Poincare- 
invariance of these potentials. 
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